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surface span, ft 

drag coefficient, Drag fore e/qS 

induced drag coefficient, Induced drag force /(/S' 

aerodynamic lift coefficient, L/qS 

induced lift coefficient, Induced lift force/§5 

zero-lift pitching-moment coefficient, Zero-lift pitching moment /(/Sc 
thrust coefficient, T/qS 

effective thrust coefficient, CY cos 6 V + »/TRCr(l — cos 6 V ) 
thrust-loss coefficient, htlCt(1 ~ cos 6 V ) 
surface chord, ft 
induced drag 

influence term as defined in equations (3.1.8) and (4.1.10) 

efficiency factor 

Hamiltonian (see eq. (2.5.1)) 

incidence angle, rad 

induced lift parameter 

aerodynamic lift force, lb 

horizontal separation between effectors, ft 

Mach number 

number of surfaces 

effective drag penalty 

dynamic pressure, lb/ft 2 

surface area, ft 2 

thrust force, lb 

free-stream velocity 

tail volume coefficient 

total lift coefficient of configuration, Weight/^S 
vertical separation between effectors, ft 
angle of attack of body reference line, rad 
thrust deflection angle, rad or deg 
average downwash angle at lifting surface, rad 
fraction of thrust recovered 

angle between gravitational horizontal and body reference line, rad 

Lagrange multiplier 

fraction of thrust loss, 1 — t/^r 
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a Prandtl coefficient 

P) normalized term, Quantity /Reference quantity 

Subscripts: 

eg center of gravity 

j surface j 

k surface k 

jk interaction terms of surfaces j and k 

v jet nozzle 

jv interaction terms of surface j and jet nozzle (see eq. (4.1.10)) 

1 wing 

2 aft horizontal tail 

3 canard 
Abbreviations: 

LOTS linear optimum trim solution 

s.m. static margin 

VLM vortex- lattice method 



Summary 

Airplane designs are currently being proposed with a multitude of lifting and control devices. Because of 
the redundancy in ways to generate moments and forces, there are a variety of strategies for trimming such 
airplanes. A linear optimum trim solution (LOTS) is derived using a Lagrange formulation. LOTS enables 
the rapid calculation of the longitudinal load distribution resulting in the minimum trim drag in level, steady 
state flight for airplanes with a mixture of three or more aerodynamic surfaces and propulsive control effectors. 
Comparisons of the trim drags obtained using LOTS, a direct, constrained optimization method, and several 
ad hoc methods are presented for vortex-lattice representations of a three-surface airplane and a two-surface 
airplane with thrust vectoring. These comparisons show that LOTS accurately predicts the results obtained 
from the nonlinear optimization and that the optimum methods result in trim-drag reductions up to 80 percent 
compared with the ad hoc methods. 

1. Introduction 

In the interest of increased economy, maneuverability, and safety, airplane designers are proposing 
configurations with more than two lifting surfaces or longitudinal-control effectors. Examples of these 
configurations are airplanes with three lifting surfaces (e.g., wing, tail, and canard) or two lifting surfaces 
and thrust vectoring. Such airplanes are capable of generating more forces than are necessary to maintain 
longitudinal equilibrium and are thus statically indeterminate. The redundant forces and moments present 
in these configurations may add many desirable capabilities to the airplanes such as attitude and flight path 
decoupling, increased damage tolerance, and enhanced aerodynamic efficiency. Some of these benefits have 
been studied by Kendall (1984), Butler (1983), Rokhsaz and Selberg (1985), and Capone and Reubush (1983). 
With the increased flexibility that these redundant forces provide, the designer faces new challenges not present 
in conventional designs. One such challenge is that the individual surface lift coefficients or thrust-vector angles 
are no longer uniquely determined by the conditions of trimmed flight. For example, the wing, horizontal tail, 
and canard of a three-surface airplane produce three independent lift forces (within aerodynamic limits). To 
maintain level trimmed flight, two constraints must be satisfied; the lift must equal the weight of the airplane, 
and the sum of the moments acting on the airplane must be zero. Thus, the number of independent forces 
exceeds the number of constraints, and hence there is an unlimited number of load distributions that can 
satisfy the trim constraints. The designer can capitalize on this indeterminateness by selecting the longitudinal 
load distribution that will satisfy the constraints and maximize or minimize a given performance function. For 
cruising flight, a logical parameter to minimize is trim drag. 

During the design process, one approach for obtaining the load distribution yielding the minimum trim 
drag for a specific flight condition is to employ nonlinear programming techniques. These numerical methods 
use multiple iterations to converge on an optimum solution. Typically, a nonlinear optimization routine 
implemented on a computer passes a load distribution to a cost-function evaluation routine. This evaluation 
routine returns the value of a cost function, which is usually based on a weighted sum of the deviation from 
trim and the trim drag as determined for the load distribution. The optimization routine compares the value 
of the cost function for the current load distribution with values from previous iterations and, based on this 
comparison, modifies the distribution in an effort to reduce the cost function. This process continues until the 
optimization routine converges on the load distribution that yields the lowest obtainable cost-function value. 
Even though this method generally yields the optimum distribution, it has several shortcomings that may 
make it unsuitable for a number of situations. For a given configuration, selecting the proper weighting terms 
in the cost function may require a fair amount of experimentation. In addition, such algorithms frequently 
need operator supervision to ensure that convergence is occurring properly. Finally, the iterative nature of this 
type of solution does not by itself provide a designer with insight into the physics governing the airplane. To 
obtain this insight, a potentially lengthy series of parametric design studies should be conducted on a given 
configuration. 

An analytical solution to the optimum trim problem that addresses many of the limitations of the nonlinear 
optimization approach has been developed and is detailed in this report. The algorithm is referred to as the 
linear optimum trim solution, or LOTS. LOTS uses a Lagrange formulation to obtain the longitudinal load 
distribution for a specific flight condition that will minimize a quadratic model of trim drag while satisfying 
linearized trim constraints. The LOTS algorithm is closed form and computationally efficient, thus making 
it compatible with a designer’s intentions and needs during the conceptual and preliminary design processes. 
The speed of the method allows timely optimization of the load distribution, whereas the closed-form nature 



permits analysis of the influence effects between control effectors. In addition to its applications in airplane 
design, LOTS may also have applications in wind-tunnel testing and in-flight operations. These potential uses 
are discussed briefly in this paper. 

The accuracy of LOTS is evaluated by comparing the trim drags obtained from its utilization with those 
obtained through a nonlinear optimization method. LOTS is also* compared with procedures that have been 
proposed as ad hoc methods of trimming redundantly controlled airplanes during preliminary design. These 
comparisons were conducted using vortex-lattice models of a general aviation class, three-surface airplane and 
a two-surface, thrust-vectoring, high-performance airplane. 

2. Derivation of LOTS 

The goal is to develop a closed-form analytical procedure for optimally scheduling the forces generated 
by the longitudinal-control effectors of a redundantly controlled airplane with respect to trim drag in steady, 
unaccelerated flight. Equations are presented that model the prime aerodynamic, propulsive, and gravitational 
effects; and methods for computing each of the relevant aerodynamic and interference terms are either described 
or referenced. The model includes vertical force and pitching- moment trim equations along with the classical 
induced drag expression. An approximate form for the propulsive thrust loss due to vectoring of thrust 
is developed for incorporation into an effective drag function. This drag function combines the effects of 
aerodynamic drag and the thrust loss due to vectoring into a single term by treating the thrust loss like an 
equivalent increase in aerodynamic drag. Influence coefficients modeling the aerodynamic interaction of lifting 
surfaces and thrust vectoring are included in the equations. The equations presented in this paper represent 
an airplane with any number of aerodynamic surfaces and a single jet-exhaust nozzle. The equations can easily 
be generalized to include the effects of multiple nozzles, but only at the expense of an increase in the number 
of terms. 

After developing the equations modeling the airplane, it is possible to solve the problem with nonlinear 
programming techniques using nonlinear constrained optimization. This approach is desirable if the equations 
have been augmented with experimental data and there is no time limitation for the computationally intensive 
calculations. An alternative is presented that is consistent with an airplane designer’s goals at preliminary 
design and with the predicted accuracy of the data available at this stage of the design cycle. In this method, 
the trim equations are linearized and a Lagrange formulation is used to minimize the drag function while 
satisfying the trim constraints. The mathematical formulation is described in this section and is applied to 
two specific airplane configurations in sections 3 and 4. 

Once an optimum longitudinal distribution of lift along the vehicle has been determined, the surface 
deflections and thrust deflection angle can be scheduled to generate the specified lift distribution. Equations 
for estimating the required stabilizer incidence and/or control surface deflections for a given required lift are 
presented in appendix A. 

2.1. Vertical-Trim Equation 

In figure 1 the aerodynamic and propulsive force coefficients acting on a multiple-lifting-surface airplane 
equipped with thrust vectoring in level, steady state flight are shown. The aerodynamic coefficients acting on 
the fuselage are combined with those of the primary lifting surface (i.e., the wing) to form a wing/body pair 
that will simply be referred to as the coefficients for the wing. To maintain level, trimmed flight it is necessary 
for the airplane to maintain constant altitude. This condition will be satisfied if the vertical aerodynamic and 
propulsive forces balance the gravitational force acting on the airplane. This balance is achieved for an airplane 
with n aerodynamic surfaces if 

n 

[(Cij 4- kjC T 6 v ) cos(a - ej - 9) + Cpj sin (a - ej - 0)] Sj + C t sin(<5 v + 9) - W = 0 (2.1.1) 

j=l 

In this and in the following equations, a caret over a variable (e.g., S) indicates that this variable has been 
nondimensionalized by a reference quantity. In this paper, the planform area of the wing is used as the reference 
area and the mean aerodynamic chord of that surface is used as the reference length. The term W is the total 
lift coefficient^ of the configuration; for the purposes of this paper in which only equilibrium cruising flight is 
considered, W is approximated by Weight /qS. The effective lift coefficient of each surface is expressed as the 
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sum of the lift coefficient of the surface in the absence of thrust- induced effects (Cij) and the change in the 
lift of the surface due to supercirculation (kjCx6 v ). This change due to supercirculation is a result of the jet 
wake acting to modify the circulation about the surface; this effect can be quite significant for surfaces located 
in close proximity to a jet wake. More details on supercirculation and its effects can be found in appendix B. 
The terms a, e, and 0 are the angle of attack of the body reference line, the average downwash angle at the 
lifting surface, and the angle between the gravitational horizontal and the body reference line, respectively. 
The thrust coefficient Cy is given as T/qS . 


C L,3 + k 3 S v c T 



(a) Nondimensionalized forces acting on a multisurface, thrust- vectoring airplane. Fuselage contribution is included in terms for 
main lifting surface. In equation (2.1.1), the term €* (not shown here) represents the change in local angle of attack due to 
downwash. 



(b) Relevant lengths. All lengths are measured in a body fixed axis system with origin located at 
aerodynamic center of wing. 

Figure 1. Relevant forces and lengths. 

The formulation of LOTS requires that the trim equations be linear. Equation (2.1.1) can be linearized 
when considering cruising flight by making several standard assumptions regarding the relative magnitude 
and contributions of the drag, thrust, and angular terms. Since this linearization is somewhat configuration- 
dependent, it will be demonstrated when specific examples are presented in sections 3 and 4. 

2.2. Moment-Trim Equation 

The moment-trim equation, dictated by the condition that the pitch rate and acceleration be zero, requires 
the sum of the moments about any point to be zero. If the aerodynamic center of the wing is chosen as the 
reference point, this condition will be satisfied if the forces and moments are governed by 
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C m ,o + W (jt C g 
+ sin(a - 6j 


cos 6 + 3 cg sin — Cj (i v sin ^ + z v cos <!>*,) — ^ j (C^ j + fyCr^u) [cos(a — ej)£j 

j = 2 

)2jj + Cx>j jcos(a — 6j)?j - sin(a — ^Sj = 0 (2.2.1) 


Here, C m , 0 is the zero-lift pitching moment, £ C g and z cg are the normalized horizontal and vertical separations 
between the center of gravity and the aerodynamic center of the wing, respectively, l v and z v are the separations 
of the exhaust nozzle, and ij and Zj are the separations of the aerodynamic trim surfaces. These separations 
are measured in the body axis system as indicated in figure 1. Like equation (2.1.1), equation (2.2.1) is easily 
linearized for cruising flight using common assumptions and will be demonstrated for specific configurations. 


2.3. Induced-Drag Equation 

Durand (1935) and Butler (1983) express the induced drag of a system of three elliptically loaded lifting 
surfaces as 

^11^1 ^22^2 + ^33^3 + 2<712 LjL 2 + 2ai 2 L 1 L 3 + 2a 2 3 ^2^3 

e l&l e 2 & 2 e 3^3 eyi ^1&2 ei2 ^ 3 e 23 &2 ^3 

The subscripts 1, 2, and 3 refer to the wing and two trim surfaces, respectively. The Prandtl coefficients 
account for the flow interference between surfaces due to the presence of vortex systems. The efficiency factors 
€ij account for the effects of viscous flow. 

Equation (2.3.1) can easily be extended to analyze systems with any number of lifting surfaces. This 
extension is accomplished by considering the surfaces in pairs and forming the appropriate influence terms. 
The resulting equation is most easily expressed in this extended form as a summation. Thus, 



(2.3.1) 


*-=■££ 


°jk 


^ Pi *=i e jk h jh 


LjLjc 


(2.3.2) 


where n represents the number of surfaces. For the development of LOTS it is convenient to write the equation 
in nondimensional form by normalizing the terms by qS\. The general induced drag equation can thus be 
written as 

C D ,i = £ E + kjS v C T )(C LJl + k k 6 v C T ) (2.3.3) 

j= 1 k= 1 n 3 k 

As before, the effective lift coefficient of each surface is separated into a component generated by the free 
stream and a thrust- induced component. 


2.4. Thrust-Loss Expression 

When utilizing thrust vectoring to reduce trim drag, some means of accounting for the decrease in the 
available propulsive force needs to be included. This thrust loss arises from the diversion of propulsive thrust 
to generate a trim force. For the current study, this effect is accounted for by adding a penalty term to the 
aerodynamic drag function. This penalty takes the form of the coefficient of thrust loss. Thus, 

Ctl = MTL Qr(l “ cos^) (2.4.1) 

The Ct( 1 — cos S v ) term accounts for propulsive losses from the diversion of thrust from the nominal, 
undeflected position. The fraction of thrust-loss term //xl indicates the fraction of this diverted thrust that 
is not recovered through the process of thrust recovery to provide an effective propulsive force. The thrust- 
recovery process is a result of favorable interactions between the jet wake and the aerodynamic circulation 
about the airplane and can significantly reduce the thrust loss associated with thrust vectoring. Experimental 
data indicate that typical values of /xtl will be between 0 and 0.5 for moderate deflection angles (6 V < 20°) 
(Lowry et al. 1957). Additional information on the process of thrust recovery is given in appendix B. 
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2.5. Lagrangian Formulation 

The optimum longitudinal load distribution is defined as the solution to the vertical-trim and moment-trim 
equations yielding the minimum trim drag. As mentioned in the “Introduction,” the rigorous approach to 
determine this distribution would be to employ nonlinear programming techniques. This approach will result 
in the solution with the best performance, but, unfortunately, it is often too burdensome to obtain for many 
applications. 

If one is satisfied with the fidelity provided by the quadratic effective-drag model (eq. (2.3.3)) and linearized 
forms of the trim equations, a relatively simple, closed-form solution (LOTS) can be obtained by using a 
Lagrange formulation as presented by Bryson and Ho (1975). As will be shown in ensuing numerical examples, 
the performance degradation for using this simplified approach is negligible for the conditions of cruising flight. 

The Lagrange formulation involves forming the Hamiltonian Ff, which is the sum of the quadratic 
performance index and the linear equality constraints appended through Lagrange multipliers. For this study, 
the performance index used is the sum of the induced drag equation and thrust-loss equation, and the linear- 
equality constraints are the linearized vertical-trim and moment-trim equations. Thus, H is formed as 


H = C Di + Ctl + A 1V T + ^2,MT 


(2.5.1) 


where the subscripts VT and MT indicate vertical trim and moment trim, respectively. It should be noted 
that by properly reformulating the performance index of /f, this technique could be applied to other trim 
optimization problems such as minimizing control deflections or achieving a desired body attitude. The partial 
derivatives of H with respect to the control variables and the Lagrange multipliers are equivalent to the gradient 
of the performance index with respect to the control variables while satisfying the trim equations; thus, the 
sufficient conditions for a minimum are that the partial derivatives of H be equal to zero. For LOTS, the 
independent design variables are the surface lift coefficients and the thrust deflection angles. Thus, performing 
the differentiation and setting the derivatives equal to zero give the system of simultaneous linear equations 

9H/dC L x = 0 ' 
dH/dC L j = 0 


OH/ dC L n = 0 ► 

dH/dSv = o 
dH/d Ai = 0 
dH/d A 2 = 0 . 


(2.5.2) 


where the subscript n denotes the nth lifting surface. When solved, this system yields the linear optimum 
values of the lift coefficients and thrust-vector angle. 

For clarity, solutions for a three-lifting-surface airplane and a two-lifting-surface airplane with thrust 
vectoring will be illustrated in the following sections. The representative drag and nonlinear trim equations 
will be presented for these configurations. The trim equations will then be linearized and LOTS will be formed. 
The effectiveness of the solution will be evaluated by applying LOTS in numerical examples and comparing the 
results with the results obtained using a nonlinear optimization routine and several ad hoc methods currently 
in use. 


3. LOTS for Three-Lifting-Surface Airplane 

Configurations with three lifting surfaces may be attractive for use in general aviation and in transport 
airplanes because of the capability of yielding minimum trim drag at all center-of-gravity locations. To utilize 
this capability it is necessary to do some type of longitudinal load-distribution optimization. This section 
describes the development and validation of a simple method for performing the optimization. 

3.1. Development of Solution for Three Lifting Surfaces 

For a three-lifting-surface airplane without thrust vectoring, the general trim and drag equations 
(eqs. (2.1.1), (2.2.1), and (2.3.3)) are applied by setting n = 3 and assuming that 6 V = 0. Thus, these 
equations can be written as follows: 
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Vertical trim: 


3 

EPt j cos (a — Cj — 9) + C]jj sin(a — ej — 0)] Sj + Cj sin 9 — W = 0 (3.1.1) 

j = 1 


Moment trim: 


C m ,o + W (i cg cos 9 + % g sin 9 j - Ctz v 
3 

— ^2 {C/.j[cos(a — €j)lj + sin(a: — ej)zj] + Cdj[ cos(a — tj)zj — sin(a — ej)Cj] j Sj = 0 (3.1.2) 

j= 2 


Drag: 




(3.1.3) 


Here, the subscript 1 refers to the wing, 2 to the aft horizontal tail, and 3 to the canard. If the airplane is in 
cruising flight, it is reasonable to neglect the contributions of drag to the trim equations since the drag forces 
are generally much less than the lift forces and are usually accompanied by significantly shorter moment arms. 
The thrust coefficient terms can be neglected for the same reasons. In addition, it is reasonable to assume that 
the angles 0, e, and a are small, thus allowing sine terms to be approximated by the angle and cosine terms 
by 1. Using these approximations, the trim equations (eqs. (3.1.1) and (3.1.2)) can be written as 


Linearized vertical trim: 


Linearized moment trim: 


c l,i + Cl, 2^2 + C L3 Sz -W = 0 
C m ,o + Wtc g — C 12^2^2 ~ Ci 3 l 3 S 3 = 0 


(3.1.4) 

(3.1.5) 


The trim equations (eqs. (3.1.4) and (3.1.5)) are now linear and can be adjoined with the drag-due-to-lift 
equation (eq. (3.1.3)) to form the Hamiltonian 




&ik s l 


L L ^J± S ^ C LJ C L,k + Ai(Cl,i + c La s 2 + C L , 3 S 3 - W) 

j = 1 fc=l * e 3 k °3° k 


+ A 2 (Cm,o + we eg - Cl, 2?2§2 - C L<3 7 3 S 3 ) 


(3.1.6) 


Setting the partial derivatives of the Hamiltonian with respect to the independent variables and the Lagrange 
multipliers (C/,,1, Cl 2, C/,3 , Ai, and A2) equal to zero, and arranging the equations such that the influence- 
coefficient matrix will be symmetrical, gives 


Cl,i 1 


Treiifrf 

252(719 

7rei2&i02 

25^(7ij 

^136163 

1 

0 

-1 

0 

Cl , 2 


2 52fl2 

7rei2&i&2 

25252(792 

7TC22&2 

252 5.3(72 3 
7re235l^203 

s 2 

-S' 2?2 


0 

Cl , 3 

— 

25 3 (7i 3 
^13^1 ®3 

252^3^23 
7re 23^1 &2&3 

2535.3(733 

7re 33 6^ 

s 3 

-sfe 3 


0 

Ai 


1 

§2 

S 3 

0 

0 


W 

A2 - 


0 

-S 2^2 

-s 3 l 

0 

0 


Cm,o-we cg . 


(3.1.7) 


The structure of the influence-coefficient matrix provides an interesting insight into the nature of the 
interactions between surfaces. The influence matrix is composed of only three distinct types of terms. These 
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three types are area terms, tail-volume terms, and aerodynamic-influence terms. With this idea in mind, the 
system of equations can be rewritten in the form 


' c l, r 


■Ell El 2 E13 S\ 0 

-1 

0 

Cl, 2 


E \2 E22 E23 S 2 ^2 


0 

Cl, 3 

= 

Ei 3 E23 E33 S3 V3 


0 

Ai 


& 

O 

O 


W 

. A 2 . 


O 

O 

iff 

O 
1 


■ Cm,o ~ W ^cg - 


The area terms 5, are normalized areas. The tail-volume terms V, are standard measures of the moment- 
generating effectiveness of a trim surface. The influence-coefficent terms E,j are measures of the induced drag 
generated by a surface or the interference drag created between surfaces. It is significant to recognize that the 
method used to obtain these influence coefficients is at the discretion of the implementer; the terms can be 
obtained through engineering judgment, analytical means, or experimental data. This flexibility allows LOTS 
to accommodate the various types of data that may be available during the design and testing processes of an 
airplane. A relatively simple analytical method for estimating the influence terms is presented in appendix C. 

The following numerical example demonstrates the application of LOTS to a three-surface airplane. 


3.2. Numerical Example for Three-Lifting-Surface Airplane 

A three-view drawing of the example airplane and its specifications are shown in figure 2. The first step 
toward forming the solution is to determine the Prandtl and efficiency coefficients, a and e, respectively. This 
process is detailed in appendix C and the results are as follows: 


— = 1.00 

eu 


— = 0.203 
e 12 


— = 1.00 — = 1.00 

e 22 e 33 


— = 0.102 — = 0.144 

ei3 e 23 



Component 

Area, 
ft 2 

Span, 

ft 

Taper 

ratio 

Sweep of 
1 /4-chord, deg 

A 

2 

A 

i 

Wing 

167 

46.5 

0.40 

0 

n/a 

n/a 

Tail 

41.4 

13.7 

.75 

25 

1.8 

4.32 

Canard 

22.3 

10.6 

.75 

0 

-.62 

-6.0 

Center of gravity (with s.m. 

= 10%)-- 



0 

-0.15 


Figure 2. Three-view drawing and specifications of example three-surface airplane. 
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In this analysis, the coefficients are assumed constant over the range of flight conditions considered; this 
condition results in the solution taking a very simple form as will be seen shortly. In practice, these coefficients 
would vary somewhat as a function of the surface loadings. As these variations are apt to be small over the 
range of usable cruise lift coefficients, the effect should be minimal. The remaining area and tail- volume terms 
are easily obtained, thus allowing the system of simultaneous equations (eq. (3.1.8)) to be formed as follows: 


' c l, r 


- 0.0493 

0.00840 

0.00547 

1 

0 - 

-1 

0 

Cl, 2 


0.00840 

0.0348 

0.00348 

0.248 

-1.07 


0 

C L , 3 

= 

0.00547 

0.00348 

0.0165 

0.134 

0.802 


0 

Ai 


1 

0.248 

0.134 

0 

0 


w 

. A 2 . 


. 0 

-1.07 

0.802 

0 

0 . 


- C m ,o ~ W ^cg - 


With all the terms constant in the influence-coefficient matrix constant, the matrix is independent of 
the flight condition. This independence allows the influence-coefficient matrix to be easily inverted without 
knowledge of the flight condition, thus resulting in a simple schedule for the optimum values of Cl, i, C 12 and 
Cl , 3 with respect to induced trim drag. Thus, 


C L , 1 = 0.9661V - 0.0174(C m , o + M cg ) 

(3.2.2) 

C L , 2 = 0.08041V + 0.432 (C m ,o + M cg ) 

(3.2.3) 

Cl, 3 = 0.1071V - 0.671(C m , o + Meg) 

(3.2.4) 


Values are shown in figure 3 of Cl,u Cl , 2 and Cl , 3 obtained by applying the above control schedule for 
C m ,o = —0.10, a static margin that is 10 percent stable, and W ranging from 0 to 2.0. The fuselage inclination 
and control surface deflections required to produce these lift coefficients are shown in figure 4. These positions 
were obtained through the methods described in appendix A. It is interesting to note that the optimum load for 
the aft tail is downward. Although this result may not be intuitively obvious, it is also not unexpected; previous 
research by Sachs (1978), McGeer and Kroo (1983), Kendall (1984), and others has shown that the optimum 
load distribution may include a negatively loaded trim surface for statically stable two-surface configurations. 




Figure 3. LOTS surface schedule for example three-surface Figure 4. Fuselage inclination, tail incidence, and canard ele- 
airplane with s.m. = 10 percent stable and C m , 0 = - 0 . 10 . vator angle versus configuration lift coefficient for example 

three-surface airplane with s.m. = 10 percent stable and 
C m ,o = - 0 . 10 . 

3.3, Validation of Three-Lifting-Surface Solution 

To assess the validity of LOTS as applied to the example three-lifting-surface airplane, a comparison was 
made of the results obtained from LOTS, a nonlinear optimization routine, and two ad hoc methods suggested 
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as a simplified means of trimming a three-surface airplane. The nonlinear optimum solution represents the 
theoretical minimum trim drag and thus provides an ideal baseline by which to assess the performance of the 
other solutions. The nonlinear optimization routine used is the direct, sequential-search simplex algorithm 
explained and illustrated by Olsson (1974, 1975). It is extremely reliable and robust in terms of convergence, 
although it suffers from slow convergence in regions of the independent design variables with low gradients of 
the performance index. It can be argued that computer resources are much less expensive than the manpower 
required to supervise other more efficient nonlinear optimization algorithms that need frequent adjustments to 
ensure proper convergence. Therefore, the robustness and reliability of the simplex algorithm make it highly 
desirable to use. The first ad hoc distribution considered was to carry no load on the tail. This method is 
frequently advocated by proponents of canard airplanes who claim that since the aft tail usually contributes 
to trim by producing a down load, the trim drag will be reduced by unloading it. The second ad hoc method 
was to carry equal but opposite loads on the tail and canard. This method was suggested by Kendall (1984) as 
optimum for a three-surface airplane with all surfaces in the same plane and with trim surfaces of equal span. 

The vortex-lattice method (VLM) described by Margason and Lamar (1971) was used to evaluate the trim 
drag of the load distributions as determined by the four trim solutions. The VLM is a relatively sophisticated 
method of determining the lift and induced drag of a multiplane system in uniform, potential flow. The method 
consists of modeling the system with a lattice of “horseshoe” vortices, that is, applying boundary conditions 
to solve for the strengths of the vortices and using these strengths to determine the forces acting on the 
system. The VLM has been shown to underpredict induced drag slightly; but as long as the paneling scheme 
is unchanged, this underprediction will be in a manner consistent with exact solutions (Feifel 1976). Thus, the 
VLM is highly suitable for comparisons of the same system of aerodynamic surfaces at various loadings. 

Because LOTS trims the airplane according to linearized trim equations, there is a small amount of trim 
error present when this distribution is substituted into the nonlinear trim equations. As will be shown, this 
error is very small and can be neglected for most applications. However, for this study in which the nonlinear 
trim equations were used to obtain the comparison solutions, it was found that this trim error caused load 
distributions obtained from LOTS to sometimes have slightly less trim drag than the theoretical minimum. To 
eliminate this contradiction, it was decided to reduce the trim error of LOTS-obtained distributions to a level 
comparable to that of the other distributions — the sum of the squares of the vertical-trim and moment-trim 
errors less than 10~ 9 . This reduction was accomplished by constraining one of the surface lift coefficients to 
the value obtained from LOTS and using an iterative search technique to reduce the trim error by varying 
the two remaining surface lift coefficients; in essence, the problem was made determinate around the LOTS 
distribution. It should be realized that when using nonlinear trim conditions, it is necessary to use an iterative 
search technique to obtain the two ad hoc solutions as well. 

Figures 5 and 6 demonstrate the effects of the trim-error reduction for the three-surface airplane at a static 
margin of 10 percent stable and with C m ,o = —0.10. Figure 5 shows the specified lift coefficients before and 
after the trim-error reduction for the wing and canard. (The tail lift coefficient is held constant.) Figure 6 
shows the predicted trim drag before and after the trim-error reduction compared with the results obtained 
from the nonlinear optimization routine. These results are typical; the trim-error reduction resulted in only a 
small deviation from the LOTS distributions. The wing and canard lift coefficients usually moved much less 
than 1 percent from the original value. These very small deviations imply that for most applications, the trim 
error can be neglected. 

As mentioned above, the trim drag obtained through the use of the nonlinear optimization routine is the 
minimum achievable while satisfying the nonlinear trim conditions and thus can serve as a standard for judging 
the other solutions. Figures 7, 8, and 9 show the percent of increase in trim drag over the nonlinear optimum 
solution for the three simplified approaches at various static stability margins and zero-lift pitching moments. 
Figure 7 represents the airplane at a typical cruise condition with a 10-percent-stable static margin and at 
Cm,o = —0.10. Over the range of cruise lift coefficients (0.3 < W < 0.9), LOTS is only negligibly less efficient 
than the nonlinear solutions and is clearly superior to the two alternate procedures. In figure 8 the static 
stability margin of the airplane is reduced to neutral stability; again, LOTS performs almost as efficiently as 
the nonlinear solution and much better than the alternate procedures. In figure 9 the static stability margin is 
returned to 10 percent stable, but the zero-lift pitching moment has been increased to —0.20. The trends of the 
two previous figures are repeated, and LOTS performs very favorably. These three combinations of pitching 
moment and static margin serve to demonstrate that although the two ad hoc procedures have regions in which 
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Figure 5. Effect of trim-error reduction on surface lift 
coefficients. 



Figure 7. Percent of increase in trim drag over nonlinear 
optimum solution for example three-surface airplane 
with s.m. = 10 percent stable and C m ,o = -0.10. 



Figure 6. Effect of trim-error reduction on percent of in- 
crease in trim drag for example three-surface airplane 
with s.m. = 10 percent stable and C m , 0 = —0.10. 



Figure 8. Percent of increase in trim drag over nonlinear 
optimum solution for example three-surface airplane 
with s.m. = 0 and C m , 0 = —0.10. 



Figure 9. Percent of increase in trim drag over nonlinear 
optimum solution for example three- surface airplane 
with s.m. = 10 percent stable and C m , 0 = -0.20. 
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they may yield near-minimum trim drag, only LOTS achieves this minimum over the range of lift coefficients 
studied. Furthermore, the optimum methods (LOTS and the nonlinear optimum) can provide trim-drag savings 
up to 80 percent compared with the ad hoc procedures. 

It is difficult to make a direct comparison of the computational time required to obtain LOTS and 
the nonlinear optimum solution. These times will vary significantly depending on how the solutions are 
implemented. For the current example, one can get an idea of how great this difference can be by considering the 
number of calculations required to arrive at the optimum loading. For LOTS, all that needs to be done to obtain 
the optimum distribution is to substitute the appropriate numbers into equations (3.2.2), (3.2.3), and (3.2.4). 
By comparison, the nonlinear optimization routine used in this study requires about 400 iterations through 
the nonlinear trim equations and the VLM drag analysis routine to converge on the optimum distribution. 
As mentioned before, the direct, sequential-search simplex algorithm used in the nonlinear optimization is not 
necessarily the fastest such algorithm available, but even more efficient algorithms would still have to perform 
numerous iterations. 

4. LOTS for Two-Lifting-Surface, Thrust-Vectoring Airplane 

Although thrust vectoring is usually considered a means for increasing agility or improving short field 
performance, it can also be used for reducing drag in cruising flight. By taking advantage of the redundant 
forces and moments that thrust vectoring provides, it is possible to optimize the longitudinal load distribution 
between the aerodynamic surfaces and the thrust force such that the effective trim drag is less than the trim 
drag without thrust deflection. By treating the thrust force in a manner similar to the force generated by 
an aerodynamic surface, the linear optimum trim solution of the previous section has been used to develop a 
scheduling algorithm that exploits the benefits of thrust vectoring during cruising flight. 

4.1. Development of Solution for Two-Lifting-Surface, Thrust- Vectoring Airplane 

The first step in developing the method is to write the nonlinear trim and drag equations (eqs. (2.1.1), 
(2.2.1), and (2.3.3)) as appropriate for a two-surface airplane with thrust vectoring. This is done by setting 
n = 2 in equations (2.1.1), (2.2.1), and (2.3.3) with the following result: 

Vertical trim: 

2 _ 

y [(Clj + kjCx6 v ) cos(a - ej — 6) + Cjyj sin(a — ej - 0)] Sj + C t sin(<5u + 9) — W = 0 (4.1.1) 

j = 1 

Moment trim: 

C m ,o + W cos 6 + z C g sin 6 ^ — Cj (jt v sin 6 V + z v cos <5^ 

- {{Cl , 2 + k 2 C T f> v ) [cos (a - e 2 )l 2 + sin(a - e 2 )z 2 \ 

+ C D , 2 [cos {a - e 2 )z 2 - sin(a - € 2 )^ 2 ] } S 2 = 0 (4.1.2) 

Induced drag: 

2 2 c 

C D ,i = E E 7rbTS&( C Lj + WvC T ){C Ltk + k k 6 v C T ) (4.1.3) 

j - 1 k~ 1 3 k 

The ‘equations are linearized by making the same approximations regarding cruising flight as made in the 
three-surface solution regarding small angles and drag. However, unlike the three-surface solution, the thrust 
terms with components vertical to the body reference line are retained since there is now independent control 
over its direction and these terms are usually accompanied by long moment arms. These approximations yield 
the following vertical-trim and moment-trim equations: 

Vertical trim: 

(C LA + * 1 C T S V ) + (C L , 2 + k 2 C T 6 v )S 2 + C t (6 v +0)-W = 0 (4.1.4) 
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Moment trim: 


C m ,o ~ CjiySy — (Cl , 2 + + Wt C g = 0 (4.1.5) 

In level flight a and 6 as defined in figure 1 are equal, thus making it possible to approximate the fuselage 
inclination angle by subtracting the incidence angle between the wing and fuselage (i\) from the ratio of the 
wing lift coefficient (C^) to the linear lift-curve slope (Cl x ). That is, 


C L ,l 


^0,1 


-h 


(4.1.6) 


With this relation, the vertical-trim equation can now be written in the final linearized form as follows: 
Vertical trim: 


Cl, i 



+ C 1,282 + 6 v Cx(ki + &252 + 1 ) — Cxh — W = 0 


(4.1.7) 


The moment-trim equation is expressed in the necessary form by grouping the thrust-dependent terms: 
Moment trim: 


Cm,o C^Zy — C La l 2 S 2 fiyCx ^£y “h — 0 


(4.1.8) 


In addition to the induced drag of the aerodynamic surfaces, the cost function must include the thrust-loss 
term to account for the propulsive force lost because of vectoring for trim. Adding the thrust-loss coefficient 
to the drag-due-to-lift equation gives the effective drag penalty: 


p = £ £ ■^ L rr^k( c Lj + W vC T ){c Lyk + k k s v c T ) + mtlQKi - cos s v ) (4.1.9) 

j=ik=i we i kt> j° k 


Using equations (4.1.7), (4.1.8), and (4.1.9) to form the Hamiltonian and setting the partial derivatives of H 
with respect to Ci i, Cx t 2 , S v , Ai, and A 2 equal to zero gives a system of five simultaneous equations that can 
be solved to find the optimum values of Cxi, C 12 , and 6 V based on the linearized equations. Thus, 


’ Cl, 1 “ 


~En E12 E\ v G\ 0 

-1 

0 

Cl, 2 


E\2 E 22 &2v C2 G4 


0 

6 v Ct 

= 

E\ y E2y Eyy G% G$ 


0 

Ai 


G\ G2 G3 0 0 


W + C T h 

. a 2 . 


. 0 g 4 g 5 0 0 _ 


CtZ v — C m ,o ~~ we cs _ 


(4.1.10) 
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where 


G\ 

G 2 

g 3 

g 4 

G 5 

E\\ 

Ei2 

E\v 

E 22 

E 2 V 


1 + 


C r p 

ciz 


s 2 


fcl + ^2^2 + 1 

~S 2 i 2 

(t v + k 2 S 2 (. 2 ^j 

2Si*n 

nenbl 


2S 2 a\ 2 

wei 2 bib 2 


_2_ / crngxfci + cri 2 S 2 k 2 \ 

7t6i V en^i ei2i>2 / 

2S 2 a 22 S 2 

7r(’22^2 

252 / <^12^1 cr 22 S 2 k 2 \ 

nb 2 \ e l2&l e 22 b 2 S\ ) 

2 ( <7n5ifc^ 2(71252^1^2 0-22^2^2^2 ^ , ^TL 

7T en6f 612&1 i>2 e 2 2&2 / 


Unlike the three-surface solution (eq. (4.1.7)), the influence-coefficient matrix for this thrust-vectoring- 
equipped airplane contains the thrust coefficient. Because the thrust coefficient depends on the flight condition, 
the influence matrix is also dependent on the flight condition. This dependence makes it inconvenient to invert 
the influence matrix (as was shown for the three-surface solution) until the value of the thrust coefficient has 
been obtained or estimated. 

The following example demonstrates LOTS for a typical two-surface airplane with thrust vectoring. 

4.2. Numerical Example for Two-Lifting-Surface, Thrust-Vectoring Airplane 


A three- view drawing of the example airplane and its relevant specifications are shown in figure 10. Estimates 
of the aerodynamic efficiency and influence terms were obtained using the method described in appendix C 
and are given as follows: 


— = 1.00 — = 1.00 — = 0.116 
en e22 e i2 

The lift-curve slope of the wing was estimated using VLM analysis to be 3.46 per radian. In addition to the 
aerodynamic influence terms and lift-curve slope, estimates need to be obtained of the thrust /aerodynamic 
interaction terms ( ki and Analytical methods for estimating these terms are unknown to the authors, 

and little applicable experimental data have been published. To simplify the demonstration of LOTS, the 
values of k\ and &2 used in this example are zero. For the example airplane where the separation between the 
wing and the jet nozzle is significant and little interaction is expected, the actual value of fci is probably near 
zero. However, the tail surface of this airplane is located very near the jet nozzle and significant interaction 
may occur, thus making the actual value of &2 likely to be nonzero. To bound the effects of thrust, recovery 
values of 0, 0.5, and 1.0 were used for /i^L- Experimental data from Lowry et al. (1957) indicate that a value 
of jixL between 0.5 and 0 may be expected in practice. With all the terms in the influence matrix in hand, the 
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solution can be written as 


’ C L, 1 ' 


0.180 

0.0371 

0 

1.0 + 0.289C r 

o - 

-X 

0 

C L, 2 


0.0371 

0.0570 

0 

0.220 

-0.346 


0 

Cfky 

= 

0 

0 

utl /Ct 

1.0 

-1.96 


0 

Ai 


1.0 + 0.289C r 

0.220 

1.0 

0 

0 


w 

. A 2 . 


0 

-0.346 

-1.96C r 

0 

0 . 


- C m , 0 — W £ C g - 


(4.2.1) 


As pointed out in the previous section, the presence of Cj in the influence-coefficient matrix of equa- 
tion (4.2.1) makes the inversion of the influence matrix impractical until Ct has been estimated or measured 
for a given throttle setting and flight condition. An estimate of Cj can be obtained by using standard drag- 
estimation techniques and assuming that Ct = Cp. The values of Cj^i, Cl ? 2 , and 8 V obtained from the above 
solution are shown in figures 11 and 12 for C m ,o = —0.10, a static margin of 3 percent stable, /ztl = 0.5, and 
a Mach number ranging from 0.3 to 0.9. 



Component 

Area, 
ft 2 

Span, 

ft 

Taper 

ratio 

Sweep of 
1/4-chord, deg 

A 

z 

A 

i 

Wing 

400 

37.5 

0.35 

20 

n/a 

n/a 

Tail 

88.1 

14.7 

.46 

40 

0.05 

1.4 

Nozzle 

n/a 

n/a 

n/a 

n/a 

0 

1.8 

Center of gravity (with s.m. 

= 3%) -- 



0 

-0.056 


Figure 10. Three- view drawing and specifications of example two-surface, thrust- vectoring airplane. 


4.3. Validation of Two-Lifting-Surface, Thrust- Vectoring Solution 

To assess the validity of LOTS as applied to the example two-lifting-surface, thrust-vectoring airplane, a 
study similar to the one used to validate the three-surface solution was conducted. An investigation was made 
of the results obtained from LOTS, a nonlinear optimization routine, and trimming the airplane without using 
the thrust-vectoring capability (6 V = 0). As in the previous section, the nonlinear optimization routine used 
is the direct, sequential-search simplex algorithm described by Olsson (1974, 1975), and the VLM was used 
to evaluate the effective trim drag of the various load distributions. It should be noted that the VLM routine 
used in this comparison was incapable of estimating the effects of supercirculation. Consequently, any changes 
in aerodynamic lift due to thrust vectoring were not modeled and thrust recovery was included simply by 
selecting, a priori, a percentage of thrust that was recovered. A trim-error-reduction routine, identical to the 
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one described in the previous section, was included in the implementation of the LOTS solution. As before, 
this trim-error reduction had only a small effect on the LOTS distributions. 



Mach number Mach number 


Figure 11. LOTS surface schedule for example two-surface, Figure 12. LOTS surface schedule for example two-surface, 
thrust-vectoring airplane with = 0.5 for wing and thrust-vectoring airplane with /i-pL — 0-5 for jet. 

tail. 

As was done in the previous section, the distributions obtained from the nonlinear optimization routine, 
having the lowest achievable effective trim drag, were used as a baseline to evaluate the other techniques. 
Figures 13, 14, and 15 show the percent of increase in effective trim drag compared with the nonlinear optimum 
solution for LOTS when trimming while holding 6 V to zero with /ixl equal to 0.0, 0.5, and 1.0, respectively. As 
can be seen from the figures, LOTS yields nearly the same results as the nonlinear optimum solution over the 
range of flight conditions and /xxl studied. For a Mach number M from approximately 0.4 to 0.5, neither of 
the optimum methods results in a significant trim-drag reduction compared with trimming without using the 
thrust vectoring. This absence of savings is because the optimum thrust deflection angle in this Mach number 
range is nearly zero, and hence little reduction is possible. 

It is interesting to note the sensitivity of the effective trim drag to the amount of thrust recovery. When 
no thrust is lost (//xl = 0), the optimum solutions predict trim-drag reductions of over 70 percent compared 
with holding 6 V equal to zero for the typical cruise Mach number region (0.8 to 0.9). However, when half the 
deflected thrust component is lost (//XL = 0*5), the reduction in effective trim drag drops to around 5 percent. 
When none of the deflected thrust is recovered (^xl = 1*0), the reduction is only 2 percent. It appears that 
to use thrust vectoring as a viable means of reducing trim drag, the design should recover a high percentage 
of the deflected thrust. 

5. Other Potential Applications of LOTS 

The simplicity, speed, and closed-form nature of LOTS also make it suited for applications in airplane 
testing and operations. 

5.1. Experimental Testing 

Frequently, the procedure used for testing models in a wind tunnel is to fix the control effectors while 
sweeping through a range of angle of attack or sideslip. For conventional airplanes (i.e., two lifting surfaces and 
no thrust vectoring), this procedure will generally result in the collection of data from which the characteristics 
of the trimmed configuration can be easily and accurately interpolated. However, because of the redundancy 
of forces and moments, this technique may not result in such data for airplanes with multiple lifting surfaces 
or propulsive control effectors. In order to collect data that reflect the optimum trimmed aerodynamics of the 
model, it is necessary to set the surfaces and effectors in the optimum positions for a given test condition. 
LOTS provides a simple means of calculating the optimum lift coefficients of the surfaces, and the required 
surface positions can be estimated using the methods provided in appendix B. Ideally, this entire scheduling 
process would be automated on a programmable calculator or personal computer, thus placing little additional 
burden on the experimenter. 
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Figure 13. Percent of increase in effective trim drag 
over nonlinear optimum solution for example two- 
surface, thrust- vectoring airplane with /txl = 0. 



Figure 14. Percent of increase in effective trim drag 
over nonlinear optimum solution for example two- 
surface, thrust-vectoring airplane with ^xl = 0.5. 



Figure 15. Percent of increase in effective trim drag 
over nonlinear optimum solution for example two- 
surface, thrust- vectoring airplane with /xtl — 1-0. 


5.2. In-Flight Operations 

The speed with which the minimum trim-drag lift distribution can be found using LOTS could make in-flight 
trim optimization a possibility. Once the solution has been set up for a given configuration, all that is needed 
to determine the optimum distribution is the current flight condition. These data could be obtained from the 
air data system, the throttle position, and an estimate of the gross weight and center-of-gravity location. The 
solution of the set of simultaneous equations in real time should be well within the capability of currently 
available flight computers. The entire process could be implemented as part of the control system, thus making 
the optimization process transparent to the flight crew. 

6. Concluding Remarks 

The research described in this paper investigates a closed-form trim solution yielding minimum trim drag. 
The solution known as the linear optimum trim solution (LOTS) provides designers a technique to trim, 
with minimum trim drag, redundantly controlled airplanes in the conceptual and preliminary design stages 
with rather straightforward analytical formulas. Potential applications in experimental testing and in-flight 
operations are also briefly discussed. 
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The benefits of LOTS and the optimum longitudinal load distribution to produce minimum drag in general 
are investigated by trimming vortex-lattice representations of a three-surface airplane and a two-surface, thrust- 
vectoring airplane using LOTS, a nonlinear optimization routine, and ad hoc techniques. This investigation 
shows that the performance of LOTS is insignificantly less than the nonlinear optimum over the range of 
conditions studied. Furthermore, the optimum trim solutions have been shown to provide trim-drag reductions 
up to 80 percent compared with the ad hoc techniques for the example three-surface airplane and 70 percent 
for the thrust-vectoring airplane. These trim-drag savings are on the order of 0.5 to 1.0 percent of the total 
airplane drag. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
March 1, 1989 


17 


Appendix A 

Surface Deflections 

After the optimum longitudinal load distribution has been determined, it may be necessary to estimate 
the control surface deflections and incidence angles required to generate the specified lifts. There are many 
methods available to accomplish this task and one is free to select a method consistent with his objectives. 
The method presented in this appendix was chosen for its relative simplicity. Its use should be restricted to 
lift coefficients where no significant flow separation is expected. 

The lift coefficient of an aerodynamic surface (Datcom 1978) can be expressed as 

C L = C La (a + e + i-Aa 0 L) (Al) 


where 


Ci a lift-curve slope 

i incidence angle between zero-lift line and body reference line 

a angle of attack of body reference line 

A a$L change in i due to deflection of simple hinged control flap 
t local downwash angle 

By estimating the appropriate dependent variables and then choosing the proper value for the independent 
variable, equation (Al) can be used for scheduling most conventional surfaces. For instance, if the incidence 
of the main wing is fixed and if the optimum position of any camber-changing devices on the wing has been 
determined, a will be the independent variable used to achieve the desired wing lift coefficient. For trim 
surfaces with variable incidence, i will be the independent variable used to obtain the desired coefficient; and 
for surfaces with a hinged control flap, Aao L will be used. Procedures for estimating the dependent terms are 
as follows: 

The term Ci can be approximated by using the following relation from Nicolai (1975): 


Cx a = 


“In A 


2 + 


4 + 




1/2 


where 

A aspect ratio 

A t / C sweep of maximum thickness line 

The average downwash angle e at a surface can be found using the following methods (Datcom 1978). The 
term € 2 , which represents the downwash at the aft tail due to the wing (the canard contribution is negligible), 
is given as 

dei 

62 = fa ' “* 1 
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where 


^ = -4.44 \k a Ki^K h (cos A^ 4 ) 1/2 j 


K 1 1 

* A ~ A 1+A™ 


A'tR = 


10 - 3TR 


K H = 


_ i-lii 


/it vertical separation of wing and aft tail 

it horizontal separation of wing and aft tail 

TR taper ratio of wing 

(*eff = a + e 4 - i — Aao-bi 

A !/ 4 sweep of 1 /4-chord line 

The method above can be adapted to estimate the downwash at the wing ej by treating the wing as the aft 
tail, the canard as the wing, and calculating de^/da as 


^ = -4.44 \k a K tk K h (cos A 1/4 ) 1/2 ] 


where S c /S w is the canard area divided by the wing area. 

To estimate the upwash produced at the canard by the wing, the following method by Perkins and Hage 
(1949) can be used. The effect of the wing on the canard €3 is given as 

de2 

C3 = d^ aeSA 


where 


^ = 0 . 9 ( 0 . 722 )^ 2 ) 
da 


£ w is the horizontal separation of the canard and wing, and c is the length of the mean aerodynamic chord. 
Because of the coupling of a e ff i and a e ff ? 3 , several iterations of the scheduling process may be required to 
converge on a satisfactory solution. 

The term Aao^ can be approximated by using the following relations from Datcom (1978): 

Aa 0 L = -C la -±-6 f K' 


where 


C ( C \ ^ 

C[ 6 lift effectiveness of plain trailing-edge flap, C\ 4 - + C 3 (yf ) 

Cf/C = Flap chord/Total wing chord 


= 1.242 - 


0.59911 


= 12.98 + 
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C 3 

K' 

t/c 

*f 

and 

where 


= -10.53 -£LM2Zi 

flap-effectiveness factor at large deflections 
= Wing thickness/ Wing chord 
flap deflection, rad 


K f — K\ + K 2 8 f + K% 6 ^ + K 4 

Ki = 1.011 + 0.1740 ^ 

c 

Cf 

K 2 = 0.002053 - 0.02069 -f 

c 

K 3 = 0.0004845 + 0.00002479 -£■ 

o 

K 4 = 5.688 x 10 -6 - 1.217 x 10 -7 ^ 

c 


I 
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Appendix B 

Analytical Approximations for Effects of Supercirculation 

Experimental studies have shown that the deflected jet of a thrust- vectoring airplane can significantly alter 
the flow field about the airplane. The name given to this alteration is supercirculation. Supercirculation has 
two primary effects — induced lift and thrust recovery — as clarified in the subsequent discussion. 

Induced Lift 

Induced lift is a change in the lift produced by the aerodynamic surfaces of an airplane due to the presence 
of thrust vectoring. The deflected jet acts in a manner similar to a jet flap, that is, altering the circulation 
about the surfaces and thus altering the amount of lift produced. The induced lift coefficient Cip was modeled 
as a function of thrust coefficient Cp and thrust deflection angle 6 V . Figure B1 shows thrust deflection angle 
versus normalized C x (Lowry et al. 1957). This figure indicates that Cx ? p can be modeled as 

Cl,V = kS v C T 

where k is a constant depending on surface and nozzle geometry factors. 

Thrust Recovery 

Thrust recovery is the fraction of deflected thrust that is eventually recovered to provide an effective thrust 
component. The term “thrust recovery” is somewhat misleading; the process that actually takes place is a 
reduction in induced drag. This reduction is due primarily to an upwash field created in front of the surfaces 
of the airplane by the deflected jet. Figure B2 (Capone 1974) shows the effect of thrust deflection on the drag 
characteristics of an airfoil. The figure shows that increasing the thrust deflection angle decreases the drag 
coefficient. With the process of thrust recovery in mind, the model for effective thrust coefficient becomes 


Cr,e = Cr cos 8 v 4- *?tr Qt( 1 — cos <5v) 
where r ?tr is the fraction of thrust recovered. 




Figure B2. Effect of nozzle deflection on drag characteristics. 
AT = 0.7; C T = 0.125. 


Since in this analysis we are interested in the lost thrust and not in the effective thrust, we form the relation 
for the thrust-loss coefficient as 


Ctl ~Ct~ Cx, e = Ct( 1 - cos 8y - ??tr + ??tr C0S M “ PtlCt(1 ” cos8 v ) 

where /ip l = 1 — VTR and is the fraction of thrust loss. Experimental data presented by Lowry et al. (1957) 
indicate that between 50 to 100 percent of the deflected thrust can be recovered at low-to-moderate thrust 
deflection angles (< 20°); hence, /xpL will generally have a value between 0 and 0.5. 
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Appendix C 

A Technique for Estimating e and a 

Before the linear optimum trim solution can be implemented, estimates of the efficiency factors e and the 
Prandtl coefficients a must be obtained. The methods used to obtain these coefficients are independent of the 
linear optimum solution — the implementer is free to use any applicable method, be it analytical or experimental. 
The following equations provide a simple analytical means for obtaining these coefficients. 

The ejj terms can be estimated using figure Cl to obtain the vortex-induced drag factor 8 for surface i 
(Torenbeek 1982) and the relation e = 1/(1 + 8). 

According to Butler (1983), the e ^ terms can be written as 


and the a ^ terms can be estimated as 

_ b s f 
ajk ~ b w Y 


2 9_ 

byj 


1/2 1/2 
e jk = e ii X e jj 


1 + 


(£) 


- 1 / 2 ' 


bs 

bw 


°' 75 (tr~K 


where 


b s smaller of two spans 

b w larger of two spans 

g vertical gap between surfaces 

Inspection of the above relation reveals that oj and <jjj = 1. 

Rather than using the above equations to determine the efficiency factors and Prandtl coefficients of the 
example airplane in sections 3 and 4, the VLM described in section 3 was utilized. By using the VLM 
to determine both the influence coefficients and the trim drag of the configurations, a potential source of 
disagreement was eliminated. The results from the VLM were compared with the results obtained from the 
equations above and were found to be nearly identical. The application of the VLM to find the influence 
coefficients of the three-surface airplane is as follows. By analyzing the surfaces of the example airplane 
individually and in pairs, estimates of the aerodynamic influence terms were obtained directly. 

Prom the quadratic model of drag due to lift (eq. (3.1.3)), 


'L,jj 
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c L,jC L ,k 7re jk bjb k 
Hence, the E\\ term of the influence matrix in equation (3.1.8) can be found by 

°L, 11 

where Cd^h is the induced drag coefficient estimated by the VLM for the isolated wing generating the lift 
coefficient Ci y n- Similarly, 

e 22 = f^’ 22 x 2 x 5 2 


and 


C L, 22 


Em = x 2 x 5 3 
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The Ejk terms (where j k ) were found by taking the difference in induced drag between the surfaces tested 
in pairs and those tested individually and using the relation 

r _ &CD,i,jk 

3t ~cIJc^ 

The influence coefficients were determined in an identical manner for the example thrust- vectoring airplane. 
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Figure Cl. The vortex-induced drag factor <5 taken from page 494 of Torenbeek (1982). 
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